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. Abstract: In the framework of the semi-hard (/^-factorization) approach, we ana- 

lyze the various charm production processes in the kinematic region covered by the 
HERA experiments. 

1 Introduction 

At the energies of modern lepton-hadron and hadron-hadron colliders, the interaction dynamics 
is governed by the properties of parton distributions in the small x region. This domain is 
characterized by the double inequality s ^> /i 2 3> A 2 , which shows that the typical parton 
interaction scale /i is much higher than the QCD parameter A, but is much lower than the 
total c.m.s. energy y'i. The situation is therefore classified as "semi-hard". In such a case, the 
perturbative QCD expansions in a s may contain large coefficients O [\og(s/n 2 )} = O [log(l/x)] 
which compensate the smallness of the coupling constant a s (/j, 2 /A 2 ). The resummation [ [l|. ^ 
of the terms [log(l/x) a s ] n results in the so called unintegrated parton distribution T{x, k 2 ), 
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which determines the probability to find a parton carrying the longitudinal momentum fraction 



x and transverse momentum k t . If the terms [log(/z 2 /A 2 



a. 



and [log(/i 2 /A 2 ) log(l/x)a s 



are also resummed, then the unintegrated parton distribution depends also on the probing 
scale n, and will be labeled as A(x, k^,/j, 2 ). That generalizes the factorization of the hadronic 
matrix elements beyond the collinear approximation (hereafter this generalized factorization 
will be referred to as ll k t factorization" [ @, §). The unintegrated parton distributions obey 
certain evolution equations (e.g., BFKL [ [|, |5], g or CCFM [ [7|, ^ P, PH) and are related 
to the conventional DGLAP [ [TT], [L^, [13], [L4| densities once the k t dependence is integrated 



out. Nowadays, the significance of the k t factorization (semi-hard) approach becomes more 
and more commonly recognized. Its applications to a variety of photo-, lepto- and hadro- 
production processes are widely discussed in the literature [ 15], 16], [17, [18], 19], [2(| [21, 

m, mi 



Remarkable agreement is found between the data 
"Of and electro-production 



221 of D* 



and the theoretical calculations regarding photo- 
mesons, and of forward jets [ as well as for specific kinematic correlations observed in 
photoproduction of D* mesons associated with jets [ at HERA. Also in hadro-production 
0, Xc [ and J/ip [ gg, P, P at the TEVATRON good agreement 



of beauty [ 
is observed 



Xc [ El and J/iP [§§[_ 
However, for a consistent application of /c r factorization in different models, the 
unintegrated gluon distribution has to be determined in the same framework. Also the various 
approximations needed to describe the experimental data have to be carefully investigated. 
In the present paper we have attempted a systematic comparison of model predictions with 
experimental data regarding the heavy flavor production processes at HERA. 



2 The fct-factorization approach applied to charm pro- 
duction 

The production of open-flavored cc pairs in ep-collisions is described in terms of the photon- 
gluon fusion mechanism. A generalization of the usual parton model to the /cf-factorization 
approach implies two essential steps. These are the introduction of unintegrated gluon dis- 
tributions and the modification of the gluon spin density matrix in the parton-level matrix 
elements. 

Here we consider only 7*g* — > cc. Let fc 7 , k g , p c and p 5 be the four-momenta of the initial 
state photon, the initial state gluon, the final state quark and anti-quark respectively, and e 7 
and e g are the corresponding polarization vectors. The photon-gluon fusion matrix elements 
for the production of an open-flavored cc pair then reads (with a charm mass m c ): 

M{ ig ^cc) = u( Pc )^ + ju( P5 ) (1) 

The matrix-element squared for open heavy quark production has already been calculated in [ 
0, which we label CE-CCH in the following. In [ [33|] (labeled as SZ) the calculation of the 



matrix elements for open heavy quark production has been repeated. In [ [20], (labeled as 
BZ) the method of orthogonal amplitudes [ 0] was applied. When calculating the spin average 
of the matrix element squared, BZ uses L^ u for the photon polarization matrix: 



IT = = 4na [8j#£ - 4(p e A; 7 )<H /{klf (2) 
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where p e is the four momentum of the incoming electron. The expression also includes the 
photon propagator factor and photon-lepton coupling. In the calculation of CE-CCH and SZ 
the photon is treated in a similar way as was the gluon in [ H : 



GT 



e 9 e g 



K gKg/\ktg\ 



(3) 



This formula converges to the usual ^2 e M e* I/ = — g^ u when kt — > 0. 

In BZ the complete set of matrix elements have been tested for gauge invariance by substitut- 
ing the gluon momenta with their polarization vectors showing explicitly the gauge invariance 
of the matrix element in order 0(a s ). 

The hard scattering cross section for a boson gluon fusion process is written as a convo- 
lution of the partonic cross section a(x g , k t ; 7*g* — > qq) with the k t dependent (unintegrated) 
gluon density A(x,k 2 ,fi 2 ) (here and in the following k t (&t 7 ) is a shorthand notation for \k t \ 
(\kt 7 |) with k t (kt 7 ) being the two-dimensional vector of the transverse momentum of the gluon 
(photon)): 



(4) 



a = / dk t dx g A(x g} k t ,n )a(x g ,k t ; 7 V -> qq) 



with the off-shell matrix elements either given by CE-CCH, SZ or by BZ. The multidimensional 
integrations can be performed by means of Monte-Carlo technique either by using VEGAS [ 
p5| 1 for the pure parton level calculations, or by using the full Monte Carlo event generator 
Cascade [ |23|, |24|, |36|. Since it is difficult to compare the different matrix elements with each 





CE-CCH 


SZ 


BZ 


k t 7 = kt = GeV 


1 


1 


0.9 


k tl = GeV, kt = 10 GeV 


1 


1 


0.86 


fc* 7 = 10 GeV, kt = GeV 


1 


1 


0.96 


k tl = 10 GeV, kt = 10 GeV 


1 


1 


0.93 



Table 1: Comparison of different calculation of the matrix elements for 7* g* — > cc. Shown are 
the matrix elements normalized to the matrix element CE-CCH for p t c = 5 GeV, r] c = r] 5 = 
in the ep c.m.s. with a/s = 300 GeV. For the comparisons the momenta of the incoming and 
outgoing partons have been modified to satisfy the small x requirement: fc 7 = Xjp e + kt 7 and 
k g = XgP p + kt resulting in k"i 



-k 2 tl and k 2 g 



-kl 



other analytically and to prove that they agree, we have performed several numerical checks. 
In Tab. [I] we show a numerical comparison of the different matrix elements for j*g* — > cc. 
For the comparison we have chosen the transverse momentum of the charm quark to be pt c = 
5 GeV, its rapidity to be r] c = i] 5 = in the ep c.m.s. and calculated the other kinematic 
quantities accordingly. The transverse momenta of the incoming partons are as indicated in 
the table. In addition we have removed all a s dependencies from the matrix elements. In 
the calculation of the off-shell matrix elements, approximations are necessary to satisfy the 
^-factorization theorem: the gluon polarization tensor as given in eq.([3|) is applied (which 
is different to the full polarization tensor of eq.(@)), and the transverse momentum must be 
dominating the virtuality: k 2 = —kf. The last condition is essentially the small x (or high 
energy) approximation. These criteria are satisfied in the calculations of CE-CCH and SZ, 
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whereas the calculation of BZ is done using the full polarization tensor IJ lv for the photon and 
only G^ v for the gluon without applying the small x approximation to the gluon four-vector. 
To study the effect of the different approximations, we have compared the matrix elements at 
large c.m.s energies of \/s = 30000 GeV and observe good agreement for all cases. This means, 
that in the asymptotic limit of very high energies, the small x approximation and the use of G^" 
are justified. In Tab. [I] we compare the matrix elements at a c.m.s energy of ^/s = 300 GeV, 
typical for HERA experiments, for different values of the transverse photon momentum k tl 
and the incoming gluon k t . Whereas the calculations of CE-CCH and SZ agree perfectly in all 
cases, a systematic difference of the order of ~ 10% to the calculation of BZ is observed. Since 
we obtained agreement in the high energy limit, this can be attributed to the effect of the small 
x approximation. In addition BZ treat the polarization of the photon and the gluon differently, 
and therefore this effect can be quantified by observing a difference of ~ 10% between the 
rows two and three in Tab. [|. The above investigations indicate that the effects of the small x 
approximation applied at HERA energies are of the order of ~ 10%. 

It is also interesting to consider the limit k t — > of the matrix elements. To do that we 
define a reduced cross section a: 

a(k t ) = J rfLips \ME\ 2 (5) 

where we integrate over the Lorentz-invariant-phase-space (Lips) of the final state quarks. The 
matrix element \ME\ 2 is taken from CE-CCH, where we have set 16n 2 a em a s e 2 = 1. In Fig. p] 




Figure 1: The reduced cross section a(k t ) as a function of the transverse momentum k t of the 
incoming gluon for different values of the transverse momentum of the incoming photon kt 7 
(m = 1.5 GeV in (a), m = 5 GeV in (b), ^/s = 30000 GeV and a fixed x 7 = x g = 0.01 ). 



we show a(k t ) as a function of the transverse momentum of the incoming gluon k t for quark 
masses of m = 1.5 GeV in Fig. |TJa and for m = 5 GeV in Fig. [ij? using ^/s = 30000 GeV and 
a fixed = 0.01. In both smooth behavior for k t — > is observed. It is also 
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interesting to note that in all cases the cross section starts to decrease at k 2 > 4m 2 . The region 
k 2 > 4m 2 is still contributing to the total cross section significantly, showing one of the main 
differences to the usual collinear approximation, where this region is completely ignored. 



3 The unintegrated gluon distributions 

Cross section calculations require an explicit representation of the k t dependent (unintegrated) 
gluon density A(x, k 2 ,fj, 2 ). We have used three different representations, one (JB) coming from 
a leading-order perturbative solution of the BFKL equations [ [37|, the second set (JS) derived 
from a numerical solution of the CCFM equation [ 23, 24] and the third (KMR) from solution 



of a combination of the BFKL and DGLAP equations [ pq] , 

JB The unintegrated gluon density A(x, k 2 ,fi 2 ), in the approach of [ [37|, is calculated 

as a convolution of the ordinary gluon density xG(x,fi 2 ) (here we use GRV [ |3"D| ) 
with universal weight factors: 

A(x,k 2 ,fi 2 ) = f g{z,kl^)-G(-,fi 2 )dz, (6) 

J x Z Z 



G(z,k 2 ^ 2 ) = ^U2^aM{l/z)\^ 2 /k 2 )) 1 fc 2 </i 2 , (7) 



G(z,k 2 ^ 2 ) = ^I (2Ja s \n(l/z)Hk 2 / l x 2 )), k 2 > /i 2 , (8) 

where Jo and Io stand for Bessel functions (of real and imaginary arguments, re- 
spectively), and a s = 3a s /n is connected to the pomeron intercept ct(0) = 1 + A, 
with A = a s 41og2 in LO. An expression for A in NLO is given in [ [4"0| : A = 
a s 41og2 — A^a 2 . In our calculations presented here we use the solution of the LO 
BFKL equation and treat A as free parameter varying between 0.166 < A < 0.53 
with a central value of A = 0.35. 
JS The CCFM evolution equations have been solved numerically in [ E3L |2~4] using a 



Monte Carlo method. According to the CCFM evolution equation, the emission of 
partons during the initial cascade is only allowed in an angular-ordered region of 
phase space. The maximum allowed angle H for any gluon emission sets the scale // 2 
for the evolution and is defined by the hard scattering quark box, which connects 
the exchanged gluon to the virtual photon. 

The free parameters of the starting gluon distribution were fitted to the structure 
function Fz(x, Q 2 ) in the range x < 10~ 2 and Q 2 > 5 GeV 2 as described in [ 24 . 



KMR In KMR [ |38[ the dependence of the unintegrated gluon distribution on the two 
scales k 2 and /i 2 was investigated: the scale fi 2 plays a dual role, it acts as the 
factorization scale and also controls the angular ordering of the partons emitted in 
the evolution. This results in a form similar to the differential form of the CCFM 
equation, however the splitting function P(z) is taken from the single scale evolution 
of the unified DGLAP-BFKL expression discussed in [ HT|. The unintegrated gluon 



5 



density xA(x, k 2 , /i 2 ) covering the whole range in k 2 has been evaluated by [42]. 



giving: 



1-2 

to 



if k t < k 



to 



xG{x,k1 ) 

- — if kt > k t0 



(9) 



with xG(x, k 2 ) being the integrated MRST [ |43f gluon density function and f(x, k 2 , fi 2 ) 
being the unintegrated gluon density of [ ^] starting from k 2 > k 2 = 1 GeV 2 . The 
unintegrated gluon density xA(x, k 2 , fj, 2 ) therefore is normalized to the MRST func- 
tion when integrated up to k 2 . 

In Fig. ^|we show a comparison of the gluon density distributions at n 2 = 100 GeV 2 obtained 
from JB, JS and KMR as a function of x for different values of k 2 and as a function of k 2 for 
different values of x. 

From Fig. we see that all three unintegrated gluon distributions show a significantly 
different behavior as a function of x but even more as a function of k t . It will be interesting to 
see, how this different behavior is reflected in the prediction of cross sections for experimentally 
observable quantities like charm production at HERA. 



4 Numerical results and discussion 

A comparison between model predictions and data in principle has to be made on hadron level 
and only if it turns out that hadronization effects are small will a comparison to parton level 
predictions make sense. However, a full simulation even of the partonic final state, including 
the initial and final state QCD cascade needs a full Monte Carlo event generator. Such a Monte 
Carlo generator based on /c r factorization and using explicitly off-shell matrix elements for the 
hard scattering process convoluted with /c r unintegrated gluon densities is presently only offered 
by the Cascade [ |3|, 



], |24], |36[ program which uses the CCFM unintegrated gluon distribution. 
This is because only the CCFM evolution equation gives a description on how to explicitly 
build the initial state gluon radiation by applying angular ordering. Other sets of unintegrated 
gluon distribution can only be used to calculate quantities at the matrix element level, which 
can be compared to data only if the effect of hadronization and of the complete initial state 
parton cascade are insignificant. 

In the following we want to systematically compare the predictions from the /^-factorization 
approach to published data on charm production at HERA. For this we use D* photo-production 
data from ZEUS [ |44| and D* production in deep inelastic scattering from both ZEUS [ |45|] and 
HI 



46|. To do this we first calculate observables using a pure parton level calculation based 



ZEUS ( 7 p) [ H 


130 < W < 280 GeV 


Q 2 < 1 GeV 2 






< 1.5 


p t > 2 GeV 


ZEUS DIS [ [H 


0.02 < y < 0.7 


1< Q 2 < 600 GeV 2 




V 


< 1.5 


1.5 <p t < 15 GeV 


HI DIS [|6| 


0.05 < y < 0.7 


1< Q 2 < 100 GeV 2 




V 


< 1.5 


1.5 <p t < 15 GeV 



Table 2: Kinematic range of the different data used for comparison 



on the matrix element calculation of BZ including the Peterson fragmentation function [ |47| for 
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Figure 2: The kt dependent (unintegrated) gluon density at fi 2 = 100 GeV 2 as a function of x 
for different values of kf (upper part) and as a function of kf for different values of x (lower 
part) as given by JS /[|| Hj|/ (solid line) , JB [ffi/ (dashed line) and KMR [\3^, g|/ (dotted 
line). 



the transition from the charm quark to the observed D* meson, where the D* meson is assumed 
to take a momentum fraction z of the charm quark, in the cc center-of-mass (c.m.s.) frame. 
Then we compare the result with a full hadron level simulation using the Monte Carlo generator 
Cascade with the matrix element of CE-CCH. Also here the Peterson fragmentation function 
is used but now with z being defined as the light-cone momentum fraction in the center-of-mass 
system of the string connecting the charm quark with its light quark partner, as implemented 
in Jetset/Pythia [ 48fl . We choose the JS unintegrated gluon for this comparison, which is 



also appropriate for a description of heavy quark production at high energies [ 28 



Next we investigate on parton level different unintegrated gluon densities. After the optimal 
choice of model parameters has been found for the JB gluon density, giving the best possible 
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agreement with data, we show a comparison to the JS and KMR unintegrated gluon density. 
We then study the sensitivity of the model predictions to the details of the unintegrated gluon 
density, the charm mass and the scale. 

We also consider the rapidity distribution of the produced D*, which is very sensitive to the 
choice of the unintegrated gluon density and the details of the c — > D* fragmentation. 

Then we investigate the x 7 distribution, which is sensitive to the details of the initial state 
cascade. We compare the predictions from a pure parton level calculation and a full event 
simulation of Cascade with the measurements. 

At the end we show, motivated by preliminary studies of ZEUS 
sensitive to the details of the heavy quark production mechanism. 

4.1 Transverse momentum distribution of D* mesons: comparison 
of parton and hadron level 

One observable which is expected to show only little sensitivity to the hadronization and to 
the full simulation of the initial and final state QCD cascades is the transverse momentum 
p t of the D* meson in photo-production and deep inelastic scattering. In Fig. [| we show 



, predictions which are 




Figure 3: The differential cross section da/dpt for D* production: (a) in photo-production 
(ZEUS [ (b) in DIS (ZEUS [ and (c) in DIS (HI [ \@). The solid line is the 

prediction from the full hadron level simulation Cascade and the dashed line shows the parton 
level calculation. In both cases the Peterson fragmentation function has been used. 



the transverse momentum distribution of D* mesons as measured by the ZEUS [ [44], |45fl and 
HI [ [5| collaborations both in photo-production and deep-inelastic scattering. The data are 
compared to the predictions of the Cascade Monte Carlo event generator on hadron level 
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including a full simulation of the partonic and hadronic final state. Also shown is the pure 
parton level calculation using the matrix element of BZ. In both cases the transition from the 
charm quark to the observed D* meson was performed by a simple Peterson fragmentation 
function (with e = 0.06 and a c — > D* branching ratio BR = 0.26). The scale ft 2 in a s {ji 2 ) 
was set to p 2 = pf + m 2 with p t being the transverse momentum in the jg c.m.s. of the final 
charm quark state assuming m c = 1.5 GeV. The JS unintegrated gluon distribution [ 24 was 
used, with the scale /i (being related to the maximum angle) p? ~ x 7 x g s. The sensitivity to 
the details of the charm fragmentation and to the full initial state gluon cascade simulation 
can be seen by comparing Cascade with the parton level calculation. We observe, that the p t 
distribution of D* mesons both in photo-production and deep inelastic scattering is in general 
well described, both with the full hadron level simulation as implemented in Cascade and also 
with the parton level calculation supplemented with the Peterson fragmentation function. We 
can thus conclude, that the p t distribution is only little sensitive to the details of the charm 
fragmentation. 

4.2 Transverse momentum distribution of D* mesons: sensitivity to 
unintegrated gluon distributions 




P.(GeV) p t (GeV) p,(GeV) 



Figure 4: The differential cross section da/dp t for D* production: (a) in photo-production 
(ZEUS (b) in DIS (ZEUS and (c) in DIS (HI [\@). The solid (dashed, dotted) 

line corresponds to using the JS (JB, KMR) unintegrated gluon density (all calculated at parton 
level). In all cases the Peterson fragmentation function has been used. 

Since the transverse momentum distribution of D* mesons is only slightly sensitive to the 
details of the full parton cascade and charm fragmentation, we can now proceed to investigate 
the sensitivity to the choice of the unintegrated gluon distribution. In Fig. ^ we show the 
prediction for da/dpt obtained from the parton level calculation as above using the JS, JB and 
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KMR unintegrated gluon distributions in comparison with the data. Although differences are 
observed in the x g and k t distributions between the different unintegrated gluon distributions 
(see Fig. 0), it is interesting to note, that very similar predictions for the D* cross sections as 
a function of the transverse momentum p t are obtained. 

In Fig. |5| we investigate in more detail the different unintegrated gluon distributions, the 
effect of varying the A parameter (see sec. (3)) in the JB distribution and of changing the 
charm-quark mass and the evolution scale fx. The comparison is performed on the parton level 
and the Peterson fragmentation function has been used to produce the D* meson. We define 
the ratio 

_ da/dp t 
~ da re f/dp t 

where da re ^ /dp t is calculated using A = 0.35, m c = 1.5 GeV and fx 2 = s/4. 




2 4 6 8 10 12 2 4 6 8 10 12 2 4 6 8 10 12 



P,(GeV) p t (GeV) p T (GeV) 

Figure 5: The ratio R = d ^rlfj^ pt as « function of p t for D* photo-production. The reference set 
corresponds to the JB set with A = 0.35, m c = 1.5 GeV and ft 2 = s/4. In (a) the solid (dashed) 
line corresponds to the JS (KMR) unintegrated gluon distribution. In (b) the parameter A of 
the JB unintegrated gluon distribution is varied (the solid (dashed) line corresponds to A = 0.53 
(0.166),). In (c) the solid line corresponds to m c = 1.3 GeV and the dashed line to the setting 
jj? = s. 



The ratio of the differential cross section as a function of p t obtained from the parton level 
calculation supplemented with the Peterson fragmentation function using the JS (solid line) 
and KMR (dashed line) unintegrated gluon distributions is shown in Fig. [5|a. The fact that the 
ratios show different behavior is directly connected to the different slopes in the k t distribution 
of the parton densities (see Fig. 0). 



10 



In Fig. |5|& we show the ratio R for different values of the A parameter of the JB unintegrated 
gluon distribution. We observe that the ratio R varies with p t in the low p f -range but seems 
to flatten off at higher p t . For A- values larger than the reference value an increase in R is 
observed, whereas A-values below the reference value result in a decreasing R. This implies 
that the p t spectra get harder with increasing A-values and softer with decreasing A values 



compared to the reference value A = 0.35 [ 20] 



In Fig. |5|c the effects of changing the mass of the charm quark and the evolution scale are 
illustrated. A decrease of the charm quark mass leads to an increase of the ratio R in the low 
p t range whereas an increase of the interaction scale from p 2 = s/4 to p 2 = s leads to an overall 
decrease of R by 20 - 25 %. 



4.3 Inclusive distribution: comparison of parton and hadron level 

In Fig. |6| we compare the measured cross section for D*-production as function of Q 2 , W, XBj 
and zd with calculations using Cascade with the JS unintegrated gluon distribution (solid 
line) and the parton level calculation supplemented with the Peterson fragmentation function 
using the JB unintegrated gluon density (dashed line). We also show the effect of changing the 
charm fragmentation (dotted line). 

Good agreement with data is observed for both Cascade with the JS unintegrated gluon 
density as well as for the parton level calculation with the JB (A=0.35) unintegrated gluon 
density in the differential cross sections as a function of log Q 2 and log x Bj [ • The differential 
cross section as a function of W, as shown in Fig. |6]&, is well described by Cascade and some- 
what less well described by JB in the peak region, although the errors of the measurement are 
fairly large. For the energy fraction zd taken by the D* meson, presented in Fig. we observe 
a slight shift of the JB distribution towards higher zjj values compared to Cascade. However, 
the Zd distribution is sensitive to the details of the D* fragmentation, which is indicated by the 
dotted line in Fig. |]g?, which represents Cascade, but using the Jetset/Pythia fragmenta- 
tion function instead of the Peterson one. In conclusion of the above comparisons, we observe a 
good description of the differential cross sections as a function of the transverse momentum p t 
of the D* meson as well as of inclusive quantities: the unintegrated gluon distributions, which 
have been considered, are reasonable for the description of the data. 



4.4 Rapidity distribution of D* mesons: comparison of parton and 
hadron level 

In photo-production and in DIS the differential cross section da/drj, where rj is the pseudo- 
rapidity of the D* meson, is sensitive to the choice of the unintegrated gluon distribution. In 
Fig. ^ we show a comparison of da/dr] in jp and in DIS at parton level supplemented with 
the Peterson fragmentation function using the JS (solid line), JB (dashed line) and the KMR 
(dotted line) unintegrated gluon distribution. Large differences in da/drj are visible, but one has 
to keep in mind that especially the r] distribution is also sensitive to the details of the c — > D* 
fragmentation, and therefore a clear distinction of the unintegrated gluon distributions based 
on this quantity alone might be questionable. In Fig. || we show da/drj, in photo-production 
and in DIS, using the JS unintegrated gluon distribution at parton level and with the full 
simulation of Cascade. We observe, that the parton level prediction including the Peterson 
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Figure 6: Differential cross sections as measured by ZEUS. Shown is a comparison of the 
calculations using Cascade with the JS (solid line) and JB (dashed line) unintegrated gluon 
density. The dotted line shows the Cascade prediction using the Jetset/Pythia charm 
fragmentation function. 



fragmentation function is not able to describe the measurements over the full range of rj. The 
effect of a full hadron level simulation is clearly visible as Cascade provides a much better 
description of the experimental data. Here the Jetset/Pythia charm fragmentation has been 
used. 

From the above it is obvious, that the r\ distribution is sensitive to the details of the fragmen- 
tation of the charm quark into the D* meson but also sensitive to the simulation of the initial 
state QCD cascade. It also shows, that an ordinary next-to-leading calculation at parton level 
cannot be expected to describe the 77 distribution, since a full simulation of the initial state 
parton evolution obviously is important for a reasonable description of heavy quark decays. 
However, in photoproduction even the full event simulation of Cascade shows differences to 
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Figure 7: The differential cross section da/dr] for D* production: (a) in photo-production 
(ZEUS cuts), (b) in DIS (ZEUS cuts) and (c) in DIS (HI cuts). Shown is a comparison of 
the calculations using the JS (solid line), JB (dashed line) and KMR (dotted line) unintegrated 
gluon density. 



the measurement at large values of 77 > 1. These large values of r] > 1 are related to large values 
of the momentum fraction x g > 0.03 of the gluons entering the hard subprocess, as shown in 
Fig. It is intuitively understandable, that the small- x approximation becomes less reliable 
at large x g -values. In addition, values of x g > 0.03 are not constrained in the determination of 



JS as described in [ [24], |28fl . It is therefore not surprising that the description is not perfect in 
this kinematic region. 

In Fig. |H] we show the cross section da/drj of deep inelastic D* production as measured 
by HI [ |46j for different ranges in zd together with the prediction of CASCADE with the 
Jetset/Pythia charm fragmentation function and the JS unintegrated gluon density (solid 
line). Also shown is the parton level calculation supplemented with the Peterson fragmentation 
function and the JB (dashed line) unintegrated gluon density. We observe that again the 
JS unintegrated gluon distribution together with a full event simulation gives a reasonable 
description of the data of this double differential cross section. 



4.5 D* and associated jet production: comparison of parton and 
hadron level 

In the BFKL and/or CCFM equations the transverse momenta of the exchanged or emitted 
partons are only restricted by kinematics. In such a scenario, the hardest pt emission can be 
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Figure 8: The differential cross section da/drj for D* for (a) photoproduction as measured by 
ZEUS [\4$, (b) DISfrom ZEUS /gf and (c) DIS from HI The histograms show the full 

hadron level simulation from CASCADE (solid line) compared to the parton level calculation 
(dashed line), both using the JS unintegrated gluon density . 




Figure 9: The distribution ofx g in D* photo-production for r] < 1 andr] > 1 forp tD * > 2 GeV. 
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Figure 10: The differential cross section da/drj for D* production in three bins of Zo as mea- 
sured by HI [ f|6|/. The data are compared to the calculations using Cascade with the Jet- 
SEt/Pythia charm fragmentation and with the JS (solid line) unintegrated gluon density. The 
dashed line shows the parton level calculation with the Peterson charm fragmentation function 
and with the JB unintegrated gluon density. 

anywhere in the gluon chain, and needs not to sit closest to the photon as required by the 
strong q 2 ordering in DGLAP. 

Photo-production of charm is an ideal testing ground for studying the underlying parton 
dynamics, since charm quarks are predominantly produced via 7 — > cc. The observation of 
any emission (jet) with p t > p1{pf) indicates a scenario, which in DGLAP is possible only in 
a full 0(a 2 s ) calculation or when charm excitation of the photon is included. However, in kt 
factorization such a scenario comes naturally, since the transverse momenta along the evolution 
chain are not k t ordered. 

The ZEUS collaboration has measured charm and associated jet production [ |44jl. In these 



measurements, the quantity of interest is the fractional photon momentum involved in the 
production of the two jets of highest E t , which is experimentally defined as 



obs _ E lt exp(-Vi) + E 2t exp(-r] 2 ) 
1 ~ 2E e y 



with E it and 77^ being the transverse energy and rapidity of the hardest jets and y being the 
fractional photon energy. In Fig. [TTJ we show a diagrammatic representation of the different 
processes involved in charm photoproduction. Thus, if the two hardest transverse momentum 
jets are produced by the cc pair, then x° BS is close to unity (Fig. |lTJa) , but if a gluon from the 
initial state cascade together with one of the c-quarks form the hardest transverse momentum 
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Figure 11: Diagrammatic representation of the processes involved in charm photoproduction. 
(a) shows the typical boson-gluon fusion diagram, the hard partons are the cc quarks, (b) shows 
a typical resolved photon (charm excitation) diagram, the hard parton are here the eg. (c) shows 
a typical diagram in kt-factorization, any two of ccg can be the hardest partons. 



jets, then x° BS < 1 (Fig. |TI]&). In a leading order calculation using the collinear approximation 
x obs < i indicates a resolved photon like process (Fig. |TT|6). Such a scenario is obtained 
naturally also in a full NLO (0(a^)) calculation, because in the three parton final state (ccg) 
any of these partons are allowed to take any kinematically accessible value (Fig. |TT|c). In the 
k t factorization approach the anomalous component of the photon (7 — > cc) is automatically 
included, since there is no restriction on the transverse momenta along the evolution chain. 

The experimentally observed x® B spectrum (Fig. |I^c) shows a tail to small values of x^ B , 
indicating that the hardest emission is indeed not always coming from the charm quarks. In 
Fig. P^ q we show a comparison of the x° BS distribution obtained on parton level with the JB 
unintegrated gluon density. Indicated is also the contribution from events where the gluon is 
the hardest, next-to-hardest and softest parton. A significant part of the cross section comes 



from events, where the gluon is the jet with the largest transverse momentum [ p5| . In Fig. |12p 
we compare the x° BS distribution obtained at parton level from the JB unintegrated gluon 
with the one from JS. In both cases a significant tail towards small x 7 values is observed, 
however the details depend on the unintegrated gluon distribution. In Fig. Il^ c we show a 



comparison of the measurement from the ZEUS collaboration [ with the prediction from the 
full event simulation of Cascade using the JS unintegrated gluon distribution and applying jet 
reconstruction and jet selection at the hadron level. We observe a reasonably good agreement, 
showing indeed a large tail towards small x^ BS values in agreement with the observation from 
data. We can conclude that the /^-factorization approach effectively simulates heavy quark 
excitation and indeed the hardest p t emission comes frequently from a gluon in the initial state 
gluon cascade. 

Another interesting quantity is the angular distribution of resolved photon like events 
(x° BS < 0.75) compared to the direct-photon like events (x° BS > 0.75) [ |4£| ^(J. An im- 
portant difference between the two (direct (Fig. |TTla), resolved (Fig. |TT1&)) scattering processes 
is that a quark (charm in this case) is the propagator of the hard scattering in the direct photon 
like events and a gluon in the dominant resolved photon like events. The angular distribution 
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Figure 12: The differential cross section dajdx^ for Q 2 < 1 GeV 2 . In (a) the parton level 
calculation is shown using the JB unintegrated gluon density (solid line). Shown are also the 
contributions where the gluon takes the largest (dash-dotted line), next-to-largest (dashed line) 
and smallest (dotted line) transverse momentum. In (b) the JB (solid line) and JS (dashed 
line) sets are used in the parton level calculation. In (c) the measurement of ZEUS [ |^|/ is 
compared to the full event simulation obtained from Cascade using the JS unintegrated gluon 
distribution. 



of dijets with a D* ± in the final state is dependent on the type of propagator (quark or gluon) 
connecting both jets (if we neglect the case, where a parton is emitted in the rapidity range 
between the two hardest jets). In the collinear approximation, the angular distribution is de- 
termined by the matrix element (7$ — > cc in the direct case or eg — > eg in the resolved photon 
case). In the /c r factorization approach the angular distribution will be determined from the 
off-shell matrix element, which covers both scattering processes. 



In Fig. 13 we show the | cos9* \ distribution, where 9* is the scattering angle of the hard jets 



to the beam axis in the dijet c.m.s. Applying the same cuts and using the same jet-algorithm 



as in [ [49], [50], we can see that the direct photon like events give a slow rise in cross-section 
with increasing | cos 9* \ . However, the cross section of the resolved photon like events rises very 
rapidly because of the t-channel gluon exchange being a combined effect of the off-shell gluon 
and the unintegrated gluon distribution. 

Thus the partons of the initial state radiation in the fc f -factorization approach give infor- 
mation on the spins of the propagators, as well as on the parton dynamics of the underlying 
sub-processes. We specifically checked the dynamics in our approach by only allowing the final 
state partons (cc) to appear in the hard scattering (i.e. turning off the simulation of the initial 
state gluon cascade in Cascade, but keeping the final state parton shower). Here we should 
expect to see only the quark exchange kind of behavior for both direct and resolved like events, 



which is also verified as seen in Fig. 13 with dotted and dashed lines overlapping around the 



default mode with x 1 > 0.75. 
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Figure 13: The distribution in |cos#*| for resolved-photon like events (x° < 0.75) and for 
direct-photon like events (x° bs > 0.75) within the kinematic range of the ZEUS measurement [ 

w 

In order to further probe the parton dynamics in this approach, we divide the entire sample 
into two parts, one where the D* proceeds along the photon direction, (i.e. t]d* < 0) and one 
where it travels along the proton direction (i.e. rjo* > 0). If the t-channel gluon is indeed 
the dominant contribution to the cross section, then the angular distribution will be peaked 
towards the direction of the incoming photon. This asymmetry should persist (more or less) 
equally for charm and anti-charm. Such an asymmetry cannot be seen in the inclusive dijet 
sample (without separation of gluon and quark jets), because by definition the distribution 
must be t u symmetric as long as we do not attempt to discriminate different kinds of jets. 

In Fig. [TJ] we show the | cos8*\ distribution for the two cases with the D* found in the 
the photon (a) and proton (b) direction, with x° BS > 0.75 and x° BS < 0.75 in the same 
phase space as in [ HUL BtJ, except with an additional cut on average pseudorapidity of the jets, 



\r]\ < 0.1, to ensure an unbiased phase space region as discussed in [ [51|. As one can see, with 



the D* in the photon direction, there is a steep rise in the cross-section for resolved photon like 
events compared to the direct photon like events. This increase in | cos 9* \ obtained through the 
initial state gluon cascade in the fc f -factorization approach can also be interpreted as "charm 
excitation" processes. On the other hand, with the D* in the proton direction, we can only see 
a mild rise in cross-section for both direct and resolved photon like events, which shows that 
the quark exchange is the dominant contribution (in the HERA kinematic range). 
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Figure 14: |cos#*| distribution with T]d* < (a) (photon direction) and t]d* > (b) (proton 
direction) within the kinematic range of the ZEUS measurement \5C 



5 Conclusion 



Three different models based on small x resummation (BFKL and CCFM formalism) and k t - 
factorization have been studied in various aspects. Two of the approaches only deliver results 
on parton level (JB and KMR) whereas the third one (JS) has been implemented into an 
event generator providing complete simulation of the initial and final state parton shower and 
hadronization. The ability of the models to reproduce experimental data has been investigated 
for charm production in the kinematic range of HERA. The aim has been to find out whether 
optimal sets of model parameters could be found leading to satisfactory description of all data, 
but also to illustrate the sensitivity of the model prediction from variations of the various 
parameters. The unintegrated gluon distributions of the three approaches exhibit different 
behaviors as a function of x but especially as a function of k t . In spite of this we find good 
agreement between the models and data on p t , logQ 2 , logW 7 and logx^- distributions, which 
indicates that these variables are not sensitive enough to differentiate between the models. The 
less good description of the zd distribution observed for two of the models might be due to 
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different treatment of the D* fragmentation. The sensitivity to the D* fragmentation is seen 
very clearly in the case of the ^-distribution which also needs a full simulation of the initial and 
final state parton emission to give a reasonable description of the data. 

From comparisons with data on photoproduction of charm leading to jets, it became evident 
that a gluon from the initial cascade frequently produces the jet of highest p t . Considering 
the uncertainties due to parton radiation and hadronization effects the models give consistent 
results. The distribution in polar angle of the hard jets, as generated by the Cascade program, 
predicts that the gluon propagator is dominant in the hard scattering of resolved photon like 
events, leading to high pt-jets initiated by a quark and a gluon. The polar angle spectrum for 
direct photon events, however, is consistent with the propagator of the hard scattering being a 
quark, resulting in hard quark-antiquark jets. 

We have shown, that the k t - factorization approach can be consistently used to describe 
measurements of charm production at HERA, which are known to be not well reproduced in the 
collinear approach. We have also shown, that in fc r factorization, resolved photon like processes 
are effectively simulated including the proper angular distributions. The k t - factorization 
approach has become now a challenging tool to understand the underlying dynamical processes 
in high energy collisions. 
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